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Abstract 

. An important issue in the dynamics of neutron star binaries is whether tidal interaction can 

' cause the individual stars to collapse into black holes during inspiral. To understand this issue 

. better, we study the dynamics of a cluster of coUisionless particles orbiting a non-rotating black 

hole, which is part of a widely separated circular binary. The companion body's electric- and 
. magnetic-type tidal fields distort the black hole and perturb the cluster, eventually causing the 

' cluster to collapse into the hole as the companion spirals in under the influence of gravitational 

. radiation reaction. We find that magnetic-type tidal forces do not significantly infiuence the 

' evolution of the cluster as a whole. However, individual orbits can be strongly affected by 

these forces. For example, some orbits are destabilized due to the addition of magnetic-type 
tidal forces. We find that the most stable orbits are close to the companion's orbital plane and 
retrograde with respect to the companion's orbit. 

O ■ 1 Introduction 

Gravitational wave detectors will soon begin operation and may see signals from sources such 
as binary neutron stars. An important issue in the dynamics of these binaries is whether tidal 
qh| interaction between the neutron stars induces each one to collapse to form a black hole during the 

late stages of inspiral (||l| and references therein). Recent numerical simulations have shown that 
bi). the stars either have a very slight tendency towards individual collapse or are completely stable 

to such collapse (see e.g. fl], |2|, |^ ) . Given the complexity of relativistic hydrodynamic simulations, 
it is useful to study simpler models in order to understand the physical underpinnings of tidally 
induced collapse (or stability against collapse) in neutron star binaries. 

One such model was considered by Shapiro and collaborators [Q, It consists of a cluster 
of coUisionless particles orbiting a non-rotating black hole, which is part of a widely separated 
circular binary. They found that as the hole's companion body spirals in slowly, its tidal field 
induces inward collapse (into the hole) of the cluster. However, they restricted attention to the 
effect of the companion body's electric-type tidal field. In the real physical situation of binary 
neutron stars, magnetic-type tidal fields may play an important role in the stability of the stars 
1^]. The purpose of this paper is to investigate the effect of the companion body's magnetic-type 
tidal field on the evolution of the particle cluster. 

2 Framework 

Consider a binary system composed of a non-rotating black hole of mass mi and a companion 
body of mass m2. Let m = mi + m2 and /x = mim2/m. Denote the binary's orbital separation by 
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b. Define the following Newtonian quantities: the relative velocity v = and the orbital 

angular velocity 

UJN = (m/63)V2. (1) 

By assumption, m <C 6, so f 1. 

Assume the binary is in a circular orbit that decays due to emission of Newtonian quadrupole 
gravitational radiation. This means the separation b evolves as [0, 

b{t) = bo{l-t/To)'/\ (2) 

where tq = (5/256)6^ (/um^)"!. 

Since the binary is widely separated, the region near the black hole can be described by a 
tidally perturbed Schwarzschild metric. Cover this region by Schwarzschild coordinates {t,r,9,(f)). 
Additional coordinate systems that will be useful are isotropic coordinates (t, x, y, z) and spherical 
isotropic coordinates {t, R,6, (p), which are related to Schwarzschild coordinates via 

r = + mi/2i?)^, = ii sin ^ cos 0, y = ii sin sin z = Rcos6. (3) 

Place a set of test particles at random positions on the sphere r = ^ b and let their initial 
velocities be in random directions tangent to the sphere. In the test particle approximation, each 
particle moves along a geodesic of the tidally perturbed Schwarzschild metric. Set the magnitude 
of each particle's initial velocity to yield a circular orbit when m2 = 0. Determine the cluster's 
time evolution by integrating the geodesic equation and assuming the particles do not collide. 



3 Estimate of magnetic- type tidal force 

Given a 3+1 split of spacetime or a local reference frame, one can separate the ten independent 
components of the Weyl tensor into an electric- and a magnetic-type tidal tensor, each having five 
independent components (see e.g. Sees. V.A.2-3 in Q). We are interested in the companion body's 



tidal fields as seen by an inertial observer in the black hole's local asymptotic rest frame [10|. This 
observer sees an electric- type tidal field whose typical components have magnitude £ ~ m2/6^, 
where "~" means "is of the order of" . If the companion body was stationary with respect to the 
black hole, the observer would not see a magnetic-type tidal field. The motion of the companion at 
velocity v with respect to the black hole induces, via a Lorentz boost with velocity v, a magnetic- 
type tidal field whose typical components have magnitude B ^ v£ as seen by the inertial observer 

The magnitude of the magnetic-type tidal force experienced by a particle moving in the local 



asymptotic rest frame is of order nipVprB ~ mpVpVrS |12|, where rup is the particle's mass, Vp is 
the particle's velocity as measured by the inertial observer, and r is the particle's orbital radius. 
This is smaller by a factor of VpV than the electric- type tidal force rriprS. 

Now consider a particle cluster at r = re- The representative particle velocity of this cluster 
is Vp = (mxjr^^l'^. Shapiro and collaborators set = 6.9mi and mi = m2, and observe cluster 
collapse when the companion body reaches b ~ 15m |5|. At that point, VpV = (mim/rcby^'^ ~ 0.1, 
so the magnetic-type tidal forces on the particles are a significant fraction of the electric ones. 
Furthermore, the two types of tidal forces typically point in different directions. Therefore, the 
magnetic-type tidal field could have an important effect on the dynamics of the cluster. In the next 
sections, we investigate this issue. 
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4 Tidal perturbation 



The leading order electric- and magnetic-type tidal deformations of a non-rotating black hole in a 
widely separated circular binary have been calculated in |11| using black hole perturbation theory. 
The tidally perturbed metric near the black hole is given in Schwarzschild coordinates {t, r, 6, (p) by 
(see Sec. Ill in 



1 _ dt" ^) dr^ + r\de' + sin2 9dcP^) 

~ ~r^) ~ ^*)^^ + sin6'cos(26l) cos((/> - iut)d())] 
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3 sin^ cos^ (0 — ut) — 1 



+ (r^ - 2ml){d9'^ + siii^ 



dt'^ + dr^ 



(4) 



where 



W = WAT ( 1 - - 



(5) 



The metric perturbation has been written in Regge- Wheeler gauge [13] and the orbital plane of the 
companion has been taken, without loss of generality, to be the black hole's equatorial plane. The 
rotation rate uj of the companion body's tidal field, as measured by an inertial observer in the black 
hole's local asymptotic rest frame, differs from the Newtonian orbital angular velocity because 



of precession effects |11|. 

The magnetic-type tidal perturbation on the black hole is embodied in the metric components 
gte and gt^ in Setting gte = ^ = gt<f> yields the metric utilized in this metric contains only 
the electric-type tidal perturbation. In what follows, we sometimes switch off the magnetic-type 
tidal perturbation by setting gte = = gtfj,, in order to isolate magnetic- type tidal effects and to 
compare with the results in |p. 

If we ignore radiation reaction, the binary's orbit remains circular and oj is constant. In this 



situation, the vector field K 



at 



+ ujSr is a Killing field of the metric (0), and K ■ u is constant 



along each particle's geodesic worldline, where u is the particle's four- velocity. This constancy 
serves as an important check on our numerics. Indeed, we first turned off radiation reaction and 
verified K ■ u was constant to sufficient accuracy along particle orbits. 

In our actual simulations, we include the effects of radiation reaction on the binary by substi- 
tuting b{t) from Eq. into (^) and (|5|). We also replace ivt in (^) by il;{t) = jQUj{t')dt', where 
uj{t) is given by (^), (|T]), and (|2|). Furthermore, we transform (^) to isotropic coordinates {t,x,y,z) 
in order to avoid difficulties with spherical coordinates. This yields 
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where R = {x^ + y'^ + z^)^!"^ . 

The particles in the cluster follow geodesies of this metric. We numerically integrate the geodesic 
equation 

(7) 



to obtain the cluster's time evolution. Here, an overdot represents a derivative with respect to 
proper time. 



5 Numerical implementation and results 

We integrate Eq. (0) using a fourth order Runge-Kutta method with adaptive stepsize control 



|14|. At f = 0, we place the companion at 6 = 60mi and put 1000 particles at Schwarzschild 
radius = 6.9mi. We also set m\ = m^- In computing the Christoffel symbols r"o-A) we ignore 
time derivatives of h(t). This is consistent with the treatment in and is justified by the wide- 
separation / slow-inspiral assumption. 

We followed the cluster's evolution until h = 15mi [t = 126068ml). To check the accuracy of 
our code, we turned off radiation reaction (i.e. set 6 to be a constant) and verified that u ■ u and 
K ■ u were conserved along particle trajectories to a fractional accuracy of 10""^ during the entire 
evolution (i.e. from t = to t = 126068mi). 

The results of the simulation are presented in Fig. ||. The cluster's average radius (r) is plotted 
against the binary's separation b for two cases: first, the magnetic-type tidal field and its distortion 
of the black hole are ignored, i.e. the metric components gu in (^) are set to zero (dashed line); 
and second, both electric- and magnetic-type tidal fields are included, i.e. the full metric (|6|) is 
used (solid line). The dashed line reproduces the corresponding result in The dynamics of the 
cluster as a whole is quite similar in the two cases, even though individual orbits can be affected 
significantly by the magnetic-type tidal perturbation (see below). A large number of particles 
plunge into the black hole aX h 30mi, causing a sudden decrease of (r) (as seen in the figure). 

It was observed in ^ that particles initially orbiting close to the equatorial plane (which is also 
the companion's orbital plane) are more stable against collapse into the black hole than those on 
initially polar orbits. Here we verify this observation in the presence of magnetic-type tidal forces. 
For each particle that plunges into the black hole by the end of our simulation, let tpiunge denote 
the time at which the particle's orbital radius is 0.51mi (2.000196mi) in isotropic (Schwarzschild) 
coordinates, and let bpiunge = ^(^piunge)- For each particle, let t denote its initial orbital inclination 
angle with respect to the equatorial plane. This quantity is defined by 

cos i = R-^{xy - yx){±^ + + i^)~^/^ . (8) 

On an unperturbed Schwarzschild background, this definition reduces to cost = L^/L, where is 
the z-component of the particle's angular momentum and L is its total angular momentum (both 
quantities are conserved along particle trajectories). 

In Fig. |2[ we plot ftpiunge versus l for each particle in the cluster. The full metric (^) is used 
to compute the cluster's evolution. Particles that have not plunged by the end of our simulation 
are represented by squares and are placed on the x-axis. Particles that plunge in the presence of 
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Figure 1: Cluster's average radius (r) as a function of binary separation b. The solid line represents 
evolution with both electric- and magnetic- type tidal perturbations present, and the dashed line 
with only electric- type present. 



the magnetic-type tidal perturbation but do not plunge in its absence are represented by triangles. 
There are no particles that plunge in the absence of the magnetic-type tidal perturbation but do 
not plunge in its presence. The rest of the particles are represented by small circles. 

A general trend in Fig. ^ is that the smaller the value of sint, the more stable the orbit. Orbits 
with 60° ^ i ^ 120° (nearly polar) are the least stable, while orbits with l Z 165° (retrograde with 
respect to the companion's orbit and nearly equatorial) are the most stable. Note that magnetic- 
type tidal forces induce collapse in some particles that were stable in the absence of these forces. 

In order to study the effect of magnetic-type tidal forces on individual orbits, we computed the 
time difference At = tp^^^gg - tpiunge' where tp^^nge (^plunge) is the plunge time with (without) the 
magnetic-type tidal perturbation on the black hole, for particles that plunge in both cases by the 
end of our simulation. This time difference is plotted as a function of the inclination angle i in 
Fig. ^. It is clear from the figure that the magnetic-type tidal force can have a significant influence 
on individual orbits. 

For 50° ^ t ^ 140°, the values of At are roughly integral multiples of 

which is the radial oscillation period of a nearly circular orbit outside r = 6mi (see also Fig. 
For i < 30°, the plunge occurs at smaller values of b (see Fig. ^) and hence the particles experience 
stronger tidal forces. As a result, the amplitudes of the particles' radial oscillations are large and the 
formula (^) for nearly circular orbits is no longer valid. In Fig. |3|, there are no points with l ^ 150° 
because no particles with l ^ 150° plunge in the absence of the magnetic-type tidal perturbation 
(see Fig. I). 

The last few radial oscillation periods before plunge for four selected particles are shown in 
Fig. ^. The vertical axis is the Schwarzschild radius of the particle. The solid line represents 
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Figure 2: Binary separation at plunge, ftpiunge; versus initial orbital inclination angle l for each 

particle. Both electric- and magnetic-type tidal perturbations are included. Squares are particles 
that do not plunge; they are placed on the x-axis. Triangles are particles that are stable in the 
absence of the magnetic-type tidal perturbation. 
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Figure 3: Difference At in plunge times with and without magnetic-type tidal forces, plotted 
against initial orbital inclination angle i. Positive values of At correspond to a delayed plunge due 
to magnetic-type tidal forces. 
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Figure 4: Radial oscillations before plunge for four selected particles: (a) i = 81.4°, (b) i = 89.9°, 
(c) L = 89.6°, and (d) l = 142°. The vertical axis is the Schwarzschild radius of the particle and 
the solid (dashed) lines represent motion with (without) magnetic-type tidal forces. 



motion in the presence of both electric- and magnetic-type tidal perturbations, and the dashed line 
in the presence of only electric. Cases (a), (b), and (c) correspond to particles with At/Tr ~ 0, -1, 
and -3 respectively. Case (d) is an example of a particle that is stable (unstable) in the absence 
(presence) of the magnetic-type tidal perturbation. 

We varied the initial cluster radius rc to investigate how this parameter affects the system's 
evolution. Our results are displayed in Fig. |^. The cluster's average radius (r) is plotted against 
binary separation b for Vc = 6.5, 6.7, 7.0, and 7.2. We put 300 particles in the cluster for each of 
these simulations. Once again, solid (dashed) lines represent motion with (without) the magnetic- 
type tidal perturbation. Cluster collapse occurs in all cases and is qualitatively similar with and 
without magnetic-type tidal forces. 

Our simulations indicate that magnetic-type tidal forces do not significantly influence the dy- 
namics of the particle cluster as a whole, though individual orbits can be strongly affected. 

Acknowledgments 

We are grateful to Kip Thorne and Lior Burko for useful comments and conversations. This 
research was supported in part by NSF grants PHY-0099568 and PHY-9796079, and NASA grants 
NAG5-10707 and NAG5-4093. 

References 

[1] G.J. Mathews and J.R. Wilson, Phys. Rev. D 61, 127304 (2000); br-qc/991104"7 . 



7 



8 
7 
6 

> 5 
V 

4 
3 

2 

60 55 50 45 40 35 30 25 

Figure 5: Cluster's average radius (r) as a function of binary separation b for Tc = 6.5, 6.7, 7.0, 
and 7.2 (with 300 particles in the cluster). The solid line represents evolution with both electric- 
and magnetic- type tidal perturbations present, and the dashed line with only electric-type present. 

[2] S. Bonazzola, E. Gourgoulhon, and J.-A. Marck, Phys. Rev. Lett. 82, 892 (1999); 

I qc/9810"07^ . 

[3] K. Uryu and Y. Eriguchi, Phys. Rev. D 61, 124023 (2000); b-qc/9908059| . 

[4] S.L. Shapiro, Phys. Rev. D 57, 908 (1998); |gr-qc/9710094i 

[5] M.D. Duez et al, Phys. Rev. D 60, 104024 (1999); br-qc/99040n . 

[6] M. Favata and K.S. Thorne (in preparation). 

[7] P.C. Peters, Phys. Rev. 136, B1224 (1964). 

[8] C.W. Misner, K.S. Thorne, and J. A. Wheeler, Gravitation (Freeman, San Francisco, 1973). 

[9] Black Holes: The Membrane Paradigm, edited by K.S. Thorne, R.H. Price, and D.A. Macdon- 
ald (Yale University Press, New Haven, 1986). 

[10] K.S. Thorne and J.B. Hartle, Phys. Rev. D 31, 1815 (1985). 

[11] K. Alvi, Phys. Rev. D 61, 124013 (2000); |gr-qc/991211^ . 

[12] W.-T. Ni and M. Zimmermann, Phys. Rev. D 17, 1473 (1978). 

[13] T. Regge and J. A. Wheeler, Phys. Rev. 108, 1063 (1957). 

[14] W.H. Press, S.A. Teukolsky, W.T. Vetterling, and B.P. Flannery, Numerical Recipes in Fortran 
77 (Cambridge University Press, New York, 1992). 




8 



